Clarification of the Bootstrap Percolation Paradox 
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We study the onset of the bootstrap percolation transition as a model of generalized dynamical 
arrest. We develop a new importance-sampling procedure in simulation, based on rare events around 
"holes", that enables us to access bootstrap lengths beyond those previously studied. By framing 
a new theory in terms of paths or processes that lead to emptying of the lattice we are able to 
develop systematic corrections to the existing theory, and compare them to simulations. Thereby, 
for the first time in the literature, it is possible to obtain credible comparisons between theory and 
simulation in the accessible density range. 



The bootstrap percolation [l|,|2j problem has attracted 
considerable interest from a variety of scientific commu- 
nities. In essence it (or its obvious variants) is a method 
to analyze the dynamics of a system of highly coupled 
units, each of which has a state that depends on those of 
its close neighbors. Such units have been considered par- 
ticles, processors, or elements of a growing population. 
Units that become active in the underlying dynamics are 
simply removed in the bootstrap. The ensemble can un- 
dergo a transition to an "arrested" state in which all dy- 
namics is quenched, and this is (after bootstrap removal 
of units) reflected in a transition from an empty to par- 
tially filled lattice. This arrest is found to be driven by 
a long length scale. The change of state of one unit be- 
comes increasingly difficult near arrest, requiring a long 
sequential string of favorable changes in its surrounding 
units. Thereby a long length (the size of the surrounding 
region to be changed) is slaved to a long (relaxational) 
time scale. 

The simplicity of bootstrap concepts means that boot- 
strap percolation plays a canonical role in the conceptual 
framework of the dynamical arrest transition, "glassifica- 
tion" ||Q,0J and in arenas as diverse as processor arrays 
0] and crack propagation amongst others. 

We study two types of "dynamics", the well-known 
"bootstrap model" itself 0, H, E3 , an d the modified 
bootstrap El • In the former, particles are removed 
if they are surrounded by c or less neighbors, and in the 
latter they are removed if any two of its vacant nearest 
neighbors are also second neighbors to each other. The 
(random bootstrap) transition occurs as a function of c, 
system size L, and initial particle density p and occurs 
when half of the prepared initial states are empty after 
removal of all movable particles. 

In this paper we introduce both new simulation al- 
gorithms and theoretical approaches that qualitatively 
change the regimes that may be explored. The theory is 
now relevant to physically accessible length-scales and, 
using only a personal computer, the simulations can be 
extended beyond the largest scales currently accessible in 
the most advanced simulations. A most interesting out- 
come, perhaps of topical interest E3> is that we are able 
to elucidate the origin of disagreements between simu- 



lation E3 and theory EE E3 for the bootstrap 
models, and show how the two can work together more 
closely in future developments. Currently our calcula- 
tions are detailed and specific to these models, but we 
consider they contain the kernel of generality required to 
signpost the path to future developments in the whole 
arena of dynamical arrest. 

The bootstrap- type problems mentioned above fall into 
two broad "universality" classes of arrest transition 0. 
The first type is a continuous or "critical" point transi- 
tion in which progressively more particles lose motion, 
leading to growing arrested domains whose typical size 
diverges with a power law. The second type of transition 
(of interest to us here) , is more reminiscent of a first-order 
transition. There, dynamically slowed domains grow near 
arrest according to an essential singularity. Mobiliza- 
tion of these domains is dependent on rare events (we 
call these "connected holes" E3) involving specific units 
that can nucleate motion on this large length. As will 
become clear, these nuclei become highly dilute near the 
transition, the typical distance between them being the 
diverging bootstrap length. 

For such transitions the following conclusions have 
been drawn by the community. For c — d (the dimension) 
it is believed that the bootstrap length £ diverges accord- 
ing to an essential singularity £ d = exp°( d )(— A/(l — p)) 
where exp ( d ) is the exponential function iterated d times 
0^|. F° r the two dimensional square lattice c = 2, the- 
oretical calculations E3 have resulted in an elegant out- 
come; essentially what are believed to be exact results, 
lim p _! 2(1 - p) log£ = A, where A = 7r 2 /9 and 7r 2 /3 for 

£ — >oc 

conventional and modified bootstrap respectively. 

On the other hand, all attempts to obtain this asymp- 
totic result by simulation have so far failed, including ex- 
tensive calculations up to L = 128, 000, leading to spec- 
ulation that it may be relevant to particle densities very 
close to unity, and consequently system sizes that are in- 
credibly large (eg. £ - 10 47 at p = 0.995 0). Such 
lengths would be far beyond what will ever be possible, 
or indeed of interest, for Physics to explore. 

Buried within this probl em, however, a more troubling 
implication emerges jj, |l2j, |l5j that simulations and the- 
ory seem relevant to such different scales that there can 
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be little useful dialogue between them. The point is also 
well made that many of the interesting arenas of appli- 
cation involve large, but finite number of units, and the 
present theory does not seem helpful there Q{. We will 
show that this need not be the case. When the problem 
is properly placed in context of general knowledge in the 
Physics community, the previous theory is useful, elegant 
and worth developing further, and the simulations were 
correct, and worthwhile. 

We will show that the modified theory and simulation 
can be brought into agreement over all reasonable length 
scales of interest to Physics. Our results are exact also 
in the true asymptotic regime |l5j but simulations can- 
not reach there and Physics will likely not be interested 
in the outcome. We identify "holes" on the lattice as 
spaces (vacancies) into which particles can move . We 
then identify these holes as either "connected" or caged 
(disconnected) according to whether the lattice can (or 
cannot) be vacated by sequentially removing particles 
beginning from that hole. The relationship to conven- 
tional (random) bootstrap simulations (described above) 
is clear; a given system size and density must contain at 
least one connected hole for it to be vacated by random 
bootstrapping processes. Thus, the bootstrap correlation 
length £ is related to the connected hole density v via 
v = l/£ 2 . The bootstrap length is therefore the average 
distance between connected holes (these representing the 
growth "nuclei" alluded to in our introductory remarks) 
that become increasingly rare near arrest. The device of 
holes allows us to focus on the key "order" parameter, 
rather than the very populous, but irrelevant particles 
and vacancies 1171. 



The details will be presented elsewhere, so here we 
present only some results for the case d = 2. In the 
simulations we begin by creating a hole with the ap- 
propriate weight, and then populate ('grow') the con- 
figuration with particles and vacancies around this site, 
checking at each stage to see if that hole is connected, or 
trapped (ie a rattler 0) by identifying cages at that 
length. Since the typical cage size grows much more 
slowly (log(l — p)/log(p)) than the bootstrap length, we 
need check only relatively small distances before the tran- 
sition probability of bootstrap to the next largest length 
approaches unity, and the hole is connected with cer- 
tainty. This approach has significant practical advan- 
tages since it permits us to sample directly only the im- 
portant rare events (connected holes) rather than prepare 
and study very large systems. Thus, the results produced 
here require only a few hours of time on a personal com- 
puter. In Figure Hand Figure |21 we show results for the 
total connected hole density in the modified and conven- 
tional bootstrap model, and these agree, where compar- 
isons are available, with the most extensive conventional 
simulations. For example, the (o) points in the upper- 
most curve of Figure represent the hole density im- 
plied by the results in HEl (system size L = 128,000), 
while the (o) points on that same curve are from our 
importance-sampling procedure discussed above. An ad- 
ditional advantage in the simulation approach is that we 
can make direct contact with theory. 

Previous theoretical calculations 0, 0, ^| approxi- 
mate the process of simultaneous removal of particles 
on increasingly large boundary contours until reaching 
one that is entirely occupied by particles, and therefore 
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FIG. 1: Modified Bootstrap Model. The (o) points represent 
the total hole density (the dotted line is a guide to the eye). 
The symmetrically growing squares (□) are compared with 
the asymptotic result exp(— 7r 2 /3(l — p)) (lower solid line) flEf . 
Also shown are results for diffusing squares (A) and small- 
asymmetry rectangles (V)- in these cases, the lines through 
the points represent our theoretical results. 
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FIG. 2: Bootstrap Model. The (o) points represent the total 
hole density (the dotted line is a guide to the eye) and can 
be compared with the results for l/£ 2 (o) for f > 100 HI20]. 
The symmetrically growing squares (□) are compared with 
the asymptotic result v — exp(— 7r 2 /9(l — p)) (lower solid 
line) |lSj|. Also shown are results for diffusing squares (A) 
and small-asymmetry rectangles (V). 
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immovable. The theory presented here is developed by 
finding the most important 'paths' or sequences by which 
particles are removed, counting the number of holes that 
would be connected using such a set of paths. Addition of 
new path types systematically improves the result, mak- 
ing comparison to theory feasible for the first time. 

Schematically, the probability of bootstrap v = is 
represented P>o — IT^ =1 (1 — p ak ) b , where b represents the 
number of sides and a the increment on the lattice. In 
the high density limit one takes the natural logarithm 
of the product, approximates the sum with an integral 
and makes the substitution y — p k . For p — > 1 this 
leads to — b/(alnp) (dy /y)\n(l — y) ~ — bn 2 /Qa{l — p) 
[l9) . In the modified bootstrap model, a = 2 and b = 4, 
one obtains — 7r 2 /3(l — p), and it is this result that has 
been proven to be asymptotically exact From Fig- 
ure Hit is clear (as has often been reported) that there is 
no agreement between the simulated results for the total 
hole density (o) and the asymptotic result, even for the 
highest densities we can simulate. It is important to ask 
how much of this deviation comes from the limited set 
of paths included in the theory and how much from the 
fact that simulations can never reach the truly asymp- 
totic limit. Thus, we calculate exactly the probability 
for removal of concentric squares of particles, a result 
that includes all corner contributions, and is valid for all 
densities. Then, 



p^ = (i-p)ii^A cs \p) 



~( CS U^ - 



(p) = 1 - V fc+1 +2p ik+1 {2 + p) 



(1) 

4j0 6fc+l + p 8k 



Using (1 — (0 2fc + 1 ) 4 as lower bounds for the coefficients 

c fc G°) an d we nn d a modified asymptotic result for the 
hole density (1 - p)~ 5 exp(-A/(l - p)) This mod- 
ified asymptotic result is almost equal to our numerical 
solution of Equation[5]and simulation results for the sym- 
metrically growing squares process in the density range 
of interest. While these are in perfect agreement with 
each other (see (□) in Figure QJ, they are still many or- 
ders of magnitude different from the full simulated con- 
nected hole density, although there is clearly a consid- 
erable improvement over the purely asymptotic result 
cxp(-^ 2 /3(l - p)). 

Nevertheless, by permitting only such restricted paths 
we have failed to correctly identify many holes as con- 
nected. Our aim now in development of the theory is 
therefore to systematically enlarge the possible paths in 
the calculation, until we approach the full simulated re- 
sult, at each stage validating the theoretical calculation 
by simulating the same restricted set of paths. By sim- 
ulation it is easy to show that indeed paths involving 
symmetric removal are the most probable, but the con- 
straint that all particles be removed from the boundary 
in a single step results in loss of many holes. If instead the 
boundary particles are simultaneously removed only on 
adjacent sides of the square, then the emptying squares 



are also permitted to "diffuse" , that is, change their cen- 
ter of gravity during the process. This leads us to identify 
many new holes that would otherwise be wrongly identi- 
fied as disconnected. One can go one stage further, and 
permit even asymmetric paths, eventually exhausting all 
paths. However, by binning the paths adopted in the sim- 
ulation of the total connected hole density we find that 
throughout the whole emptying process, beginning from 
a connected hole, the ratio of the rectangular sides never 
exceeds 1.4, by far the greatest contribution coming from 
near-square process, providing they are permitted to dif- 
fuse. We have been able to realize this approximation in 
theory also. We define intermediate states of the system 
as the squares illustrated in Figure 03 with weight P^ . 
Paths implied by removal of particles map the growing 
vacant region between only these states, larger by one 
step at each stage. If we consider only the limited set P 1 
- P 3 , such processes correspond to growing and diffusing 
squares. Inclusion of the states P 4 and P 5 permits in 
addition "fluctuations" of the square by one additional 
layer. These local intermediate states are related by the 
coupled equations, 



(2) 



where i, j's range from 1 — > n, with n = 3 in the 
diffusing squares process, and n = 5 in its extended 



small- asymmetric rectangular version. 



(p) defines 



the probability of migration from class j to i at the fcth 
step. These equations are solved subject to the initial 
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FIG. 3: Intermediate configurations arising in evaluation of 
the sum over 'paths' that may be used to empty the lattice. 
The inner square from each of these figures is assumed to have 
been emptied of particles in the previous step, and the next 
parts to be emptied lie between this and the outer bound- 
ing square. Shaded outer regions imply a complete line of 
particles, blocking further movement in that direction. Un- 
shaded outer regions have at least one vacancy in that part 
of the perimeter. The arrows on the squares indicate possi- 
ble transitions involving the growth of two adjacent boundary 
lines by one step. This process restores that local configura- 
tion to one of the intermediate states again, with the internal 
empty square one step larger in two directions. The process 
terminates when the local configuration makes a transition 
outside of the class P 1 - P 3 (diffusing squares) or P 1 - P 5 for 
small-asymmetry rectangles. Note that in P 3 we have given 
an explicit example of the extension of the boundary implied 
by the arrow. 
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conditions, = (1 — p)Su. 

The choice of these states, and transitions between 
them, is far from trivial since we must ensure an "order- 
ing" of the removal process if we wish to use a random 
measure for the particles in calculating the coefficients 
cji (p). For the sake of simplicity, we here present only 
the coefficients of the process involving diffusing squares, 

c W) = l_2^+^- 1 J 4 1 ' 2 ) = (l-p)(l-/) 

4 2 ' 2) = 1 - 2 P k + p 2k - 2 ,c^ = 2p(l - p p k -> + p k ) 
c^=p 2k ~\ cP=/(l-/- 1 ) 
cf 3 ) =p 2 (l + 2p k - 2 -4p k - l +p 2k - !i ) 

The total bootstrap probability P^>\p), may be cal- 
culated numerically for any density, limited only by the 
precision of the computer. The same processes may be 
simulated on the computer and are in each case identical 
with the theory. Results for all are given in Figure ^ and 
Figure |21 for modified and conventional bootstrap. 

The outcome is intriguing. Diffusing squares (A), im- 
proves the comparison between theory and full simula- 
tion, and small- asymmetry rectangles (V) yields results 
that may (for the first time) begin to be credibly com- 
pared to computer simulation of the full simulated hole 
density, and implicitly the bootstrap correlation length 
in the regime where simulations can be carried out. 

We may now summarize our results both in relation 
to the bootstrap problem, and in a more broad context. 
Firstly, for the bootstrap problem itself, new computa- 
tional and theoretical approaches have enabled simula- 
tions and theory to be brought into reasonable (indeed 
arbitrarily good) agreement across a wide range of den- 
sity. The fact is that both the theory and simulation 
do not adopt the very simple asymptotic form that has 
been quoted in the literature until one reaches densities 
and length scales that are beyond the natural interest 
of physics. There nevertheless remain many areas of 
physics where extended regimes of dynamic slowing, and 
near-arrest are of great importance, and these can be 
dealt with by the methods described here. Secondly, by 
properly identifying the most probable paths (incidently 
thereby respecting the symmetry of the problem), and 
developing theory as a sum over only these paths, we 
obtain a very useful approximation, and systematic cor- 
rections around it. In essence this amounts to a sort of 
"mean- field" approximation in the path-integral (sum), 
more familiar in field theory as the optimal instanton tra- 
jectory |23. Inclusion of small asymmetry is equivalent 
to the 'shape' fluctuations included in next to leading 
order in such calculations. This is more than an anal- 
ogy; the bootstrap process of this discontinuous transi- 
tion produces an essential singularity precisely because, 



underlying it, is the physics of complex activated pro- 
cesses. 

In this second point, but from a broader perspective, 
we have linked the whole bootstrap endeavor (and cru- 
cially those areas of physics for which it is considered rel- 
evant) to an arena of physics that is already somewhat 
explored, and opened the pathway to numerous develop- 
ments in the theory of dynamical arrest, many of which 
will immediately suggest themselves to the reader. 
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